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We present the observation of glasslike dynamic correlations of mobile mercury ions in the ionic 
conductor Cu2Hgl4, detected in both NMR and nonlinear conductivity experiments. The results 
show that dynamic cooperativity appears in systems seemingly unrelated to glassy and soft arrested 
materials. A simple kinetic two-component model is proposed, which seems to provide a 
description of the cooperative ionic dynamics. 
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In glass-forming materials particles increasingly move 
together as the glass transition is approached [HHf ■ Such 
cooperativity is also found in other arrested systems 0- 
Q and seems to be intimately connected to the slow dy- 
namics. Here we report on the observation of large-scale 
dynamic correlations in a distinctly nonglassy system — 
the conductive phase of the ionic conductor Cu2HgLt. 
Using carefully designed nuclear magnetic resonance ex- 
periments we prove that mercury ions are the main con- 
tributors to conduction (establishing C^HgLi as the first 
known mercury conductor), and show that mercury dif- 
fusion is anomalous. These results urge for a more de- 
tailed examination of ionic motion. Therefore nonlinear 
conductivity measurements are used as a probe for dy- 
namical heterogeneity, revealing a characteristic correla- 
tion time scale. To explain the cooperativity we pro- 
pose a simple model related to previous work on glasses 
0, HI j with two essential ingredients — disorder and ex- 
istence of two kinds of particles, slow (copper) and fast 
(mercury). We compare the results with recent studies 
of arrested and glass-forming materials 0, thus es- 
tablishing an unexpected connection between seemingly 
different fields. 

Cu2Hgl4 used in experiments was in powder form, 
synthetized according to standard procedure X-ray 
diffraction at 300 K showed no appreciable contamination 
with iodides, and all applied techniques (NMR, DSC, 
conductivity) saw a sharp transition at T c = 344.7 K, 
providing further evidence of phase purity. Free induc- 
tion decays were used to record NMR line shapes, while 
a recovery sequence was employed in the 63 Cu relaxation 
measurements. Conductivity was measured in a two- 
contact cylindrical cell with graphite electrodes, using a 
low distortion voltage source and lock-in amplifier. At all 
temperatures, sample resistance was above 1 Mil. Low 
frequency (7 Hz) conductivity agrees quantitatively with 
previously published values jlOL . The third harmonic 
current j'3 provided nonlinear conductivity; the heating 
contribution to jz was estimated to be small due to large 
sample resistance, and more importantly, uniform over 
the employed frequency range (linear conductivity is not 
peaked). Instrumental harmonic distortion effects were 
also negligible between ~ 10 Hz and ~ 20 kHz. All mea- 



surements were reproducible after several temperature 
cycles across T c . 

Before studying cooperative ionic motion in C^HgLt, 
we must identify the charge carriers and nature of the 
insulator-conductor transition at T c . Ever since the dis- 
covery of ionic conduction in Cu2Hgl4 



121, it has re- 



mained unclear which ion species predominantly carries 
current in the conducting phase above 

T c [HEl. Here 

we obtain direct proof of mercury motion from NMR 
experiments — a substantial motional narrowing of the 
mercury line in the conductive phase [Fig. HJa)]. In 
contrast, the copper line is broadened, indicating qua- 
sistatic disorder. This suggests that the transition is not 
melting of the copper sublattice (unlike the related Cul 
Q), but rather an order-disorder transition [l5[ with the 
copper ions remaining virtually static. Mercury motion 
is then enhanced in the changed energy landscape above 
T c . Slow copper diffusion and disorder are essential for 
explaining mercury dynamic cooperativity, so we perform 
NMR line shape and relaxation measurements on copper 
to provide a microscopic picture of the transition. 

The structure of Cu2Hgl4 in the ordered phase below 
T c contains 8 tetrahcdral positions per unit cell, with only 
3 occupied by copper or mercury [13J. Thus one expects 
that a relatively small activation energy is neccessary to 
create point defects by moving ions from " regular" to nor- 
mally empty tetrahedral positions. NMR enables us to 
follow these motions in an ion-specific way, and observe 
how they lead to a transition to the disordered phase. 
Naturally abundant copper nuclei are quadrupolar and 
thus sensitive to electric field gradients (EFGs) present 
in the material. This is already obvious in the NMR 
line shapes — the local environment of copper ions doesn't 
have cubic symmetry, leading to a quadrupolar splitting 
of the line. While all Cu + which are on regular sites have 
roughly the same local environment, ions on normally 
vacant sites (defects) experience much larger EFGs and 
their NMR lines are substantially shifted and broadened 
[l6j| . The signal from defects thus contributes as a low, 
broad background easily separable from the narrow line 
of "regular" Cu + . The integral of the narrow line can 
be used to obtain the relative number of defects. This 
number is expected to grow anomalously fast close to an 
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Figure 1. (color online) a) NMR lines of copper and mer- 
cury, below and above T c . Frequency is relative to Larmor 
frequencies, 145.6 MHz for 65 Cu and 91.2 MHz for 199 Hg. 
Pulse excitation width is ~ 150 kHz. Copper lines are nor- 
malized to have the same integrals, while the low-temperature 
mercury line is multiplied by 10 after normalization. Signifi- 
cant motional narrowing is observed for mercury above T c . b) 
Spin- lattice relaxation measurements for 63 Cu. Left insert is 
frequency dependence at 280 K. Lines are fits obtained from 
a superposition of phonon and two defect diffusion processes 
(mercury and copper). Temperature ranges where each of the 
processes comes into play are indicated. Right insert is order 
parameter in dependence on reduced temperature. 



order-disorder transition, as is indeed observed. One can 
even see power- law behavior close to T c , described rea- 
sonably well with a mean-field approach. We define an 
order parameter £ = (n; — n v ) / (nj + n v ), with ni and n v 
the number of regular copper atoms and copper defects, 
respectively. The transition being first order, £ drops 
abruptly to zero above T c . A Landau expansion yields 
£ — £ c ~ e* 3 close to T c , with £ c the order parameter at 
T c , e = (T c — T) /T c the reduced temperature and /3 the 
critical exponent: f3 = 1/2, comparing well with the ex- 
perimental value of 0.58 (Fig. [lb right insert). At T c the 
lattice reorganizes and the distinction between "defects" 
and "regular" copper ions is lost, as all Cu + positions 
are equally probable [13] • However, no motional narrow- 
ing of the copper line is observed above T c . Instead, the 
line broadens due to larger EFGs caused by electrostatic 
disorder. 

To learn more about defect dynamics we measure cop- 



per spin- lattice relaxation [Fig. QJb)]. Diffusing defects 
create fluctuating EFGs, which influence the relaxation. 
In addition to a conventional Raman phonon mechanism 
[H , we observe effects of both mercury and copper defect 
diffusion on the relaxation rate below T c . The contribu- 
tion from one defect type is 0, [l6| 
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with Sujq a temperature-independent quadrupolar cou- 
pling constant, Th hopping time and u>l Larmor fre- 
quency. The hopping process is thermally activated fl9| , 
with Th = Toe Eh / kT , where 1/tq is the attempt frequency 
and Eh the hopping activation energy. The numbers of 
defects also follow Arrhenius-type laws except close to 
T c . Combining the temperature and frequency depen- 
dences of the relaxation rate, with attempt frequencies 
estimated from Raman spectroscopy [2j , we obtain Eh 
for both copper and mercury defects [21( . The values are 
4900 ± 100 K and 2220 ± 50 K for copper and mercury, 
respectively. Thus already below T c mercury has a sig- 
nificantly lower hopping activation energy than copper. 
Microscopic reasons are as of yet unclear. 

Above T c the mercury diffusion rate increases for an or- 
der of magnitude and the 199 Hg line becomes motionally 
narrowed. The line shape is well fitted by a Lorentzian 
curve and the hopping time can be extracted from the 
linewidth using Aw = (Awrj) 2 Th, where Awrj is the static 
linewidth (below T c ). The simple formula is valid for 
AujQTh <C 1, so we have taken into account corrections 
for finite Th where neccessary [22I ]. Employing the Ein- 
stein relation for mobility, we can try to calculate the 
conductivity from extracted hopping times: 
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where 2e is the charge of Hg + ions, n their number 
density (~ 5 x 10 21 cm~ 3 ) and L a hopping distance 
of the order of the interatomic spacing (~ 1 A). If we 
now take this conductivity and compare it to the mea- 
sured dc values, we observe the first sign of anomalous 
behavior: in a region ~ 30 K above T c , gq is signifi- 
cantly larger than <Jdc (Fig. Hk). Thus relation ©, valid 
for simple stochastic motion of ions, doesn't correctly 
predict the long-time transport. One may ask if this is 
due to the existence of some new, intermediate timescale 
above r/j ~ 1 ^s where motional correlations arise, or just 
well-known short range correlation effects quantified with 
the Haven ratio (23|. To resolve the question, we mea- 
sure the nonlinear conductivity 173 (cj), defined by j = 
aiE + a 3 E 3 + in dependence on frequency uj (Fig. [2p). 
Although dynamical correlation effects often bear small 
influence on linear response, they are intimately related 
to nonlinear susceptibilities. A quantitative measure is 
the four-point correlation function [24H26j]. C4 (y, t) — 
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(/(x,0)/(x + r,t)/(x + y,0)/(x + y + r,t)) x (with/ 
a suitable dynamic parameter, e.g. intermediate scat- 
tering function), representing the correlation of time 
changes at different points in space. Thus if many 
ions move synchronously on a characteristic timescale 
Tcorr, C4 (y, t) will have a peak at T corr . Generalized 
fluctuation-dissipation theorems connect X4j the spatial 
integral of C4, with the corresponding nonlinear suscepti- 
bility [2i| , making dynamical correlations measurable. In 
contrast to the case of a dielectric (or magnetic) material, 
where one measures the response of dipoles to an external 
field, we detect the response of mobile charges, and the 
natural response function is 03 instead of the susceptibil- 
ity X3- A lot of activity is currently aimed at modeling 
dynamical heterogeinity in soft and glassy systems (27| . 
but experimental data are still scarce — the first report on 
nonlinear susceptibility of a glass-former (glycerol) only 
appeared recently (H. Here we see similar effects, but 
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Figure 2. (color online) a) DC conductivity (circles) is smaller 
than conductivity predicted from mercury NMR (squares), in- 
dicating anomalous behavior. The discrepancy between con- 
ductivities closely follows the number of correlated ions (di- 
amonds, estimated from nonlinear conductivity and NMR) 
b) nonlinear conductivity in dependence on frequency. Peaks 
clearly show the existence of a cooperativity time scale. Non- 
linear response below T c is negligible. Lines are guides to the 
eye. 



in a rather unexpected material. Characteristic correla- 
tion timescalcs are revealed through peaks in 03 (u>) at 
frequencies cuT corr ~ 1 , and the relative number of corre- 
lated ions, N corri can be estimated from integrals of the 
peaks [Fig. [UJa)]. In our work we focused on 03 as it 
is sensitive to correlated motion; we note, however, that 
the new correlation timescale also affects the linear con- 
ductivity <7i (w). A shouldcrlikc feature is visible at the 
frequency l/T corr in the frequency spectrum of a\\ this 
strengthens the analogy with supercooled liquids, where, 
similarly, the peak of the nonlinear response occurs close 
to the characteristic relaxation frequency visible in linear 
response. 

Except very close to the transition, r corr is substan- 
tially longer than the mercury hopping time Tu- Thus a 
simplistic conduction model can be used to explain the 
discrepancy between <7n and a^c- We assume that mer- 
cury ions move vigorously most of the time, but some- 
times get constrained to small volumes. NMR lines of 
these ions are broad and do not contribute to the prin- 
cipal narrow line. Occasionally several trapped ions ar- 
range favourably, and leave the "trap" together. Thus 
the effective number of charge carriers is diminished and 
the characteristic correlation timescale appears. This is 
essentially a "cooperatively rearranging regions" (CRR) 



28|] . A simi " 



scenario, well known in glass science 
lar mechanism was also proposed for colloidal gel relax- 
ations [H, and seems to offer a good phenomenological 
explanation of our data. Contrary to glass-forming liq- 
uids, where N corr has no effect on adc, here the corre- 
lations influence it. We obtain direct experimental evi- 
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Figure 3. (color online) Stimulated spin echo NMR measure- 
ments at 380 K, giving evidence of mercury ion trapping at 
the characteristic timescale r corr . Insert is raw measurement, 
the square denoting the zoomed-in segment where dynamic 
trapping effects are visible. Main graph is compensated for 
spin-lattice relaxation, showing only the correlation contri- 
bution. Full line is from simulation, and the dotted line a 
stretched exponential fit. 
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dence for this model from a different NMR experiment 
on mercury — stimulated spin echo (SSE) measurements 
(Fig. [3]). Moving spins experience much smaller aver- 
age local fields than trapped ones, leading to a difference 
in spin decoherence times. This can be exploited to se- 
lectively excite and detect only ions moving at a given 
moment. The SSE sequence [4j is perfectly suited for 
such an experiment (30| . After correcting for spin- lattice 
relaxation, we can directly observe how, of all ions mov- 
ing at the time of excitation, a sizeable fraction becomes 
trapped after a time A ~ T corr (Fig. [3]). This experiment 
provides an absolute scale for N corr , and we can make a 
comparison with the difference between oo and Ode [Fig. 
[2ja)]. The agreement is gratifying, both in absolute scale 
and temperature dependence, implying that the differ- 
ence can be attributed to a diminished effective number 
of carriers, confirming the phenomenological model. 

However, microscopic questions remain: what causes 
confinement, and how are correlated jumps performed? 
To answer them, we propose a very simple mechanism, 
related to investigations of spin glasses (essentialy a limit- 
ing case of the Edwards- Anderson Hamiltonian with dif- 
fusion 0, [3l|). Aside from disorder, the basic require- 
ment is the existence of two kinds of atoms in the mate- 
rial, with different diffusion coefficients, and short-range 
interactions. In C^HgLj this is realized with copper 
and mercury, on a fixed iodine background. If we take 
low-temperature activation energies to be representative, 
we can conclude that mercury diffuses 10 3 to 10 4 times 
faster than copper in the interesting temperature range 
where correlations appear. As slow copper ions move 
around, they occasionally form compartments with sev- 
eral trapped mercury ions inside. The compartments can 
then "open" due to cation rearrangement — once a path is 
open, many fast mercury ions use it sequentially to empty 
the compartment. Such behavior is indeed observed in 
a two-dimensional random walk simulation. Simulations 
were run on a square lattice with periodic boundary con- 
ditions and initially randomly placed ions. In every step 
the mercury ions moved in random (allowed) directions, 
while the copper ions moved with a certain probability 
D* (which is essentially the ratio of copper and mercury 
diffusion coefficients). In the course of simulation large 
mercury "islands" form and dissolve (Fig. 21 insert). To 
calculate the four-point correlation function Xi, we used 
the persistence function [32[, defined as n r (to) = I if 
nothing has happened on site r for t < to, and n r (to) = 
otherwise. X4 1S calculated as the variance of the autocor- 
relation of n r (t) , evaluated at all mercury sites [33| . A 
characteristic correlation timescale is revealed (Fig. @|, 
and the curves qualitatively follow the CRR prediction 
32| . The only parameters we have to set are D* and ef- 



fective particle concentration, taking care that the num- 
ber of vacant sites is above the percolation threshold [13] ■ 
For realistic concentrations and D* , correlation times be- 
come about lO 3 ^, in fair agreement with experiment. 
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Figure 4. (color online) Some results of the 2D simulation. 
Inserted frame at 2000 steps shows formed mercury islands 
(D* = 0.0022, simulation box 300x300 cells). A correlation 
timescale is nicely visible in the time dependence of the four- 
point correlation xt, f° r several D* . 



The SSE decay curve (Fig. [3J can also be predicted sur- 
prisingly well. 

The model is " minimal" , in the sense that we obtain 
dynamical heterogcinity with the minimum number of 
assumptions. This hints at a considerable universality of 
such correlations. Important effects have, however, been 
neglected: the iodine lattice potential, short-range elec- 
trostatic correlations, electronic dynamics and phonons. 
Correct temperature behavior cannot be obtained with- 
out taking them into account. In contradistinction to 
glass-forming materials, no slowing down of the corre- 
lated dynamics with decreasing temperature is observed 
in Cu2Hgl4; we suspect that this is due to the addi- 
tional periodic potential, which tends to restore an or- 
dered state (and succeeds at T c ). Also, we believe that a 
full three-dimensional simulation would show similar co- 
operative behavior (with modified vacancy numbers due 
to a lower percolation threshold), but this needs to be 
proven. More elaborate simulations are needed to better 
understand these issues. 

From our results we can conclude that ingredients 
needed for large-scale motional correlations are quite 
ubiquitous, so it is reasonable to believe that ionic co- 
opcrativity is important for many other systems as well. 
In disordered ionic conductors it might offer a more con- 
vincing explanation of nonlinear response than standard 
hopping models (35j . opening up new perspectives for 
studying ion dynamics. Even more important is the con- 
nection with arrested materials, which shows that dy- 
namical correlations are more universal than previously 
thought. The observed interplay between lattice poten- 
tial and dynamical heterogeneity is very interesting in 
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itself and could provide a unique possibility for exploring 
the emergence of glasslikc correlations. 
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SUPPLEMENTAL MATERIAL 

IONIC DIFFUSION AND COPPER 
SPIN-LATTICE RELAXATION 

Information on the diffusion of ions and pretransitional 
increase of copper defect numbers is mainly obtained 
from line shape and spin-lattice relaxation measurements 
on copper. The total relaxation rate is a sum of the 
phonon contribution and the two defect parts: 

yn = 1/Tt\ ph + l/T x \ Cu _ def + l/Tt\ Hg _ def (3) 

Quadruploar spin 3/2 copper nuclei couple to electric 
field gradients due to phonons and defects; the phonon 
part is well approximated by 1/Ti\ h = CT 2 above the 
Debye temperature (with C a constant) In the sim- 
plest approximation of point charges and exponential au- 
tocorrelation, the relaxation rate due to one kind of de- 
fects is given by Eq. (1) in the paper. Hopping times 
and defect numbers are temperature-dependent and in 
principle Arrhcnius-like, except copper defects close to 
the order-disorder transition. A linear scale plot of spin- 
lattice relaxation vs. temperature below T c (Fig. [5]) 
clearly demonstrates that two distinct defect contribu- 
tions are neccesary to reproduce the temperature depen- 
dence of the relaxation. We must therefore fit a superpo- 
sition of the form of ([3]), with many adjustable parame- 
ters. To constrain the fit, we used several additional re- 
sults. We easily determined the number of copper defects 
using measured integrals of the NMR line, as described 
in the paper. An additional field-dependent relaxation 
measurement was performed at 280 K to increase the 
reliability of the remaining fitting parameters — two acti- 
vation energies for mercury (defect formation and jump) 
and one for copper (jump), and three multiplicative con- 
stants; field and temperature dependecies were fitted 
simultaneously. Attempt frequencies were constrained 
to one order of magnitude using Raman spectroscopy 
data (where anomalous low-frequency modes are iden- 
tified with the attempt frequencies [2|), and agree with 
typical values for similar compounds (e.g. Cul [3[). 

The activation energies resulting from the fit arc given 
in the paper. The approximation of independent de- 
fects is most probably adequate for temperatures not very 
close to T c , but the reliability of our fitting model can be 
questioned near T c , as we suppose that only the number 
of copper defects grows dramatically. However, the line 
shape measurements indicate that indeed copper defects 
are much more prevalent than mercury — intensity of the 
63 Cu line drops sharply, while 199 Hg shows no appreciable 
change (within the relatively large error margin) . Above 
T c the relaxation rate is significantly higher than in the 
low temperature phase, due to the large number of fast 
mercury ions. If we assume that the relaxation model of 
Eq. (1) in the paper is at least approximately valid above 



T c , clearly the contribution of mercury ion diffusion will 
dominate the relaxation rate (as we know from the ab- 
sence of motional narrowing in Cu that copper motion is 
much slower, and wot/j 3> 1). We can thus use mercury 
hopping times obtained directly from NMR line shapes 
to make a comparison with the 63 Cu relaxation. In this 
limit l/ 63 Ti is roughly proportional to n/r^ ~ rjrj (with 
n the number of mobile ions, and <7o the conductivity 
as defined in Eq. (2) in the paper). When l/ 63 Ti, <jq 
and (Jdc are plotted on top of each other, we see that 
the relaxation rate follows the conductivities within the 
error margin (Fig. [3] insert), confirming our qualitative 
analysis. 
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Figure 5. 63 Cu spin-lattice relaxation rate below T c , on a 
linear scale. Raman phonon process and two defect diffusion 
contributions are plotted separately. Insert is relaxation rate 
above T c (circles), which is seen to roughly fall between ao 
(squares) and (Jdc (dotted line). 



STIMULATED ECHO MEASUREMENTS ON 
MERCURY 

Motional narrowing enables us to make selective NMR 
measurements on moving mercury ions. The spin-spin 
relaxation time of stationary ions, T2 j t rap , is of the order 
of 27r/Aajrj ~ 10 /is, while the fast ions have some ten 
times longer relaxation times (Fig. [SJi). If we use a stim- 
ulated spin echo sequence^ (Fig. EJd) we can exploit this 
difference. The first two pulses excite only spins which 
are moving at the moment of excitation (if T^trap "C 5 
and r5 ~ T2j as t) and thus 'tag' them magnetically. Af- 
ter a diffusion time A the third pulse is applied, causing 
a stimulated echo. However, tagged spins which have 
been trapped during the time A and are stationary at 
the moment of the application of the third pulse, have 
shorter Ti ~ Ti.tra-p and do not contribute to the echo 
signal. Thus the echo amplitude should be anomalously 
small for diffusion times comparable to and larger than 
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Figure 6. a) Spin-spin relaxation for mercury at 380 K, ob- 
tained using a conventional spin echo sequence, b) The stim- 
ulated spin echo (SSE) sequence. Radio-frequency pulses and 
spin echo signal are represented schematically. 

the dynamic correlation time r corr . The size of this 'dip' 
in diffusion time dependence provides an estimate of the 
absolute number of correlated ions. Such behavior is ob- 



served in simulation as well — even the same functional 
dependence (a stretched exponential) can be fitted to ex- 
perimental and simulation results. It is possible to obtain 
a good agreement between experiment and simulation in 
both correlation time and fraction of correlated ions (for 
a vacancy fraction of 0.4 and effective copper diffusion 
coefficient D* = 0.0008), but the simplifying assumption 
of single copper and mercury jump probabilities makes 
the stretching parameter of the simulation curve some- 
what closer to 1. This could be quickly amended by in- 
troducing a jump probability distribution 'by hand', but 
was deemed physically untransparent. The agreement 
between simulation and experiment is already quite im- 
pressive in the simple version, and more sophisticated 
models are needed to account for the finer effects. In 
absence of trapping, the echo decay should be a simple 
exponential in A, with the decay time equal to the spin- 
lattice relaxation time T\ . The spin-spin relaxation is sin- 
gle exponential (Fig. EJi), the line shape doesn't change 
noticeably for any A, and 199 Hg is a spin 1/2 nucleus. 
Hence the deviations from an exponential decay in the 
SSE experiment can only be due to dynamical trapping. 
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